Penrose Process in Kerr- Taub- NUT Spacetime 
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Abstract 

Penrose process on rotational energy extraction of 
the black hole in the Kerr- Taub- NUT spacetime is stud- 
ied. It has been shown that for the radial motion of 
particles NUT parameter slightly shifts the shape of 
the effective potential down. The dependence of the 
extracted energy from compact object on NUT param- 
eter has been found. 
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1 Introduction 

At present there is no any observational evidence for the 
existence of gravitomagnetic monopole, i.e. so-called 



NUT (Newman, Unti & Tamburino - iNewman et al 
(|l963l )) parameter or magnetic mass. Therefore study 
of the motion of the test particles and energy extraction 
mechanisms in NUT spacetime may provide new tool 
for studying new important general relativistic effects 
which are associated with nondiagonal components of 
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Morozova and Ahmedovll2009l . where solutions for elec- 
tromagnetic waves and interferometry in spacetime 
with NUT parameter have been studied.). Kerr-Taub- 
NUT spacetime with Maxwell and dilatio n fields is re- 
cently investig ated by Aliev et al. (I2008I). In our pre- 
cedin g papers ( Morozova et al.ll2008 ; Abdujabbarov et aH 
20081 ) we have studied the plasma magnetosphere 
around a rotating, magnetized neutron star and charged 
particle motion around compact objects immersed in 
external magnetic field in the presence of the NUT pa- 
rameter. 

Penrose process (see, e.g. Penrose! 19691) for the ex- 
traction of energy from rotating black hole is based 
on the existence of negative energy orbits in the er- 
gosphere, the re gion bounded by th e event horizon and 
the static limit. Gariel et al. I d2010h have examined the 
possibility that astrophysical jet collimation may arise 
from the geometry of rotating black holes and the pres- 
ence of high-energy particles resulting from a Penrose 
process rather than from effect of the magnetic fields. 
Detailed study of the energetics of the Ke rr-Newman 



black hole by the Penrose process is given bv lBhat et al 
( 19841) . Energetics of a rotating charged black hole in 5- 



di mensional supergravity has b een recently considered 
bv lPrabhu and Dadhichl |201fi ). 

The geodesies of the Kerr-Taub-NUT spacetime 
share many of the properties of the ergosphere in the 
field of a magnetic monopole. A thorough discussion 
and comparison of these orbits can be found in pa- 



pers ( Zhnmerman and Shahirll989t Kagramanova et al.^ 
2010ft . 

Here we first study motion of the test particles 
around rotating compact object with nonvanishing 
NUT parameter. The effective potential of the radial 
motion of the test particles around rotating compact 
object is investigated in the presence of the NUT pa- 
rameter using the Lagrangian formalism. Finally we 
study the energy extraction mechanism from compact 
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object through Penrose process in the Kerr-Taub-NUT 
spacetime. 

The outline of the paper is as follow. In the Sec. [2] 
we study the ergosphere and motion of the test par- 
ticles in the Kerr-Taub-NUT spacetime. The Sec. [3] 
is devoted to study the energy extraction mechanisms 
through Penrose process in the Kerr-Taub-NUT space- 
time. The concluding remarks are given in the Sec. 2] 

Throughout the paper, we use a space-like signature 
(— , +, +, +) and a system of units in which G = 1 = c 
(However, for those expressions with an astrophysical 
application we have written the speed of light explic- 
itly). Greek indices are taken to run from to 3 and 
Latin indices from 1 to 3; covariant derivatives are de- 
noted with a semi-colon and partial derivatives with a 
comma. 



2 Ergosphere around compact object in 
Kerr-Taub-NUT spacetime 

We consider electromagnetic fields of compact astro- 
physical objects in Keer-Taub-NUT spacetime which in 
a spherical coordinates (ct, r, 9, <j>) is descr i bed by the 



metr ic (see 
20031) 



Dadhich and Turakulov 2002; Bini et al 



ds 2 



A 



a 2 sin 2 i 



dt 2 



-dr 2 + Y>d9 2 



1 

2 

— [Ax - a(S + ax) sin 2 9] dtdip 
-i[(£ + a X ) 2 sin 2 0- X 2 A] dtf, 



(1) 



where parameters X, A and x arc defined by 

S = r 2 + (1+acosO) 2 , 
A = r 2 - 2Mr- I 2 + a 2 , 
X = a sin 2 9 - 2lcos9, 

M is the total mass, a is the specific angular momen- 
tum, and I is the NUT parameter of the central object. 

The spacetime (|TJ) has a horizon where the 4 velocity 
of a co-rotating observer turns to null, or the surface 
r = const becomes null: 



r+ = M + \JM 2 + P 



(2) 



The static limit is defined where the time-translation 
Killing vector £r t \ = d/dt becomes null (i.e. goo = 0) 
and static limit of the black hole can be described as 



lies outside the horizon. The region between the two is 
called the ergosphere, where timelike geodesies cannot 
remain static but can remain stationary due to corota- 
tion with the BH with the specific frame dragging an- 
gular velocity at the given location in the ergosphere. 
This is the region of spacetime where timelike particles 
with negative angular momentum relative to the BH 
can have negative energy relative to the infinity. 

In FigQ] the dependence of the shape of the ergo- 
sphere from small dimensionlcss parameter I = l/M is 
shown. From the dependence one can easily see that in 
the presence of the NUT parameter radius of the event 
horizon becomes larger. However a relative volume of 
the ergosphere is decreased. 

Due to the existence of an ergosphere around the 
black hole, it is possible to extract energy from it by 
means of the Penrose process. Inside the ergosphere, it 
is possible to have a timelike or null trajectory with neg- 
ative total energy. As a result, shoot a small particle A 
into the ergosphere from outside with energy at infinity. 
When the particle is deep down near the horizon, let it 
to explode into two parts, B and C, one of which at- 
tains negative energy relative to infinity and falls down 
the hole but other part escapes back to radial infin- 
ity by conservation of energy with energy greater than 
that of the original incident particle. This is how the 
energy could be extracted from the hole by axial accre- 
tion of particles with the suitable angular momentum 
and I parameters. Consider the equation of motion of 
such negative energy particle. The Lagrangian for this 
particle can be written as: 

2£ = - - (A - a 2 sin 2 9) i 2 + ^-r 2 + £(9 2 
£ v ; A 



[A% - a(S + ax) sin 2 9] i<p 
[(£ + ax) 2 sin 2 #-x 2 A] cp 2 , 



(4) 



and according to (j3J the generalized momenta are given 
by 



Pi 



'A - a 2 sin 2 0) i 



(A X - o(E + a X ) sin 2 9)<p = E, 



-P<p = 



= --(A X -a(Z + a X )sm 2 9)i 



--((£ + ax) 2 sin 2 (?-x 2 A)^ = £, 



= M + \lM 2 + l 2 -a 2 cos 2 9. 



(3) 



-Pr 

-Pe = 



S . 

A r ' 
= -£0, 



(5) 

(6) 
(7) 
(8) 



Considering only the outer horizon, r + and static 
limit, r st , it can be verified that the static limit always 



where superior dots denote differentiation with respect 
to an affine parameter r. (The conservation of pt and 
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Ptp follows from the independence of the Lagrangian 
on t and ip which, in turn, is a manifestation of the 
stationary and the axisymmctric character of the Kcrr- 
Taub-NUT geometry.) 

The Hamiltonian for the test particle in spacetime 
fTJ) is given by 



% = pit + p v (p + p r f + peO — C 



._(A- a W 



' 2 + I' 2 



-- [A X - o(E + a X ) sin 2 0] tip 
-^[(E + ax) 2 sin 2 0- X 2 A] p\ 



(9) 



and from the independence of the Hamiltonian on t , 
one can get 



2U = 



A 



£0 



j2 



1 

£ 



+ -(A X -a(£ + ax) sin 2 



•r 



- (A X - o(E + a X ) sin 2 



t 



1 



+ — ((£ + ax) 2 sin 2 9 - x 2 A) ip 



= Et- Lp + —r 2 + £0 2 
A 



5 = const. 
We may, without loss of generality, set 

6 = 



for time-like geodesies, 
for null geodesies. 



(10) 



(11) 



Then from the expressions and ([6|) one can easily 
obtain the expressions for ip and t as 



<P 



1 

A 



(A - a 2 sin 2 0) ( X E - L) - aY^E sin 2 9 
£ sin 2 9 



(12) 



A£L 



Ax - o(E + ax) sin 2 
| (£ + QX ) 2 sin 2 0-x 2 A 

£ (A X - a(£ + ax) sin 2 0) 
x ((A -a 2 sin 2 9)( X E - L) -aZE sin 2 0)] 



(13) 



It was first shown bv lZimmerman and Shahirl dl989h 
for spherical symmet ric case (NUT spacetime) and later 
by Bini et al. ( 20031 ) for axial symmetric case (Kerr- 
Taub-NUT spacetime) that the orbits of the test par- 
ticles are confined to a cone with the opening angle 
given by cos0 = 2EI/L. It also follows that in 
this case the equations of motion on the cone depend 



on I only via I 2 (jBini et al.l 12003c lAbduiabbarov et al 
20081 ). The main point is that the small value for the 
upper limit for gravitomagnctic moment has been ob- 
tained by comparing theoretical results with experi- 
mental data as (i) I < 10~ 24 from th e gravitational 
microlensing ( Rahvar and Habibi 20041 ). (ii) I < 1.5 ■ 
10 -18 from the interferometry exper i ment s on ultra- 
cold atoms ( Morozova and Ahmedov 2009f) . (iii) and 
similar limit has been obtained fr om the experiments 



on M ach-Zchnder interferometer (jKagramanova et al 
20081 ). Due to the smallness of the gravitomagnetic 



charge let us consider the motion in the quasi-equatorial 
plane when the motion in direction changes as = 
7r/2 + S6(t), where 59(t) is the term of first order in 
I, then it is easy to expand the trigonometric func- 
tions as sin0 = 1 - <50 2 (t)/2 + <D(59 4 {t)) and cos0 = 
60(t) - 0(59 3 (t)). Now inserting the expressions (fI2J) 
and (p~3|) to the equation (fT0|) , and neglecting the small 
terms 0(59 2 (t)) 1 one can easily obtain equation for the 
radial motion as follows: 



£f 2 



, 2 (Mr + / 2 ) 



2(L - aEy 

r l + l z 

+E 2 (r 2 +a 2 + l 2 ) - L 2 



8 A. 



(14) 



As we have noted, 6 = for null geodesies and equation 
for radial motion (IT4l) becomes 



£f 2 = 2(L-aE) : 



(Mr + f 
r 2 + 1 2 



+E 2 (r 2 +a 2 +l 2 )~L 2 . (15) 



Hereafter, it is more convenient to distinguish the 
geodesies by the impact parameter D = L/E. First, 
we study the geodesies with the impact parameter 

D = a, when L = aE, (16) 

and equations (fT2"j) . ([T3]) and (fT5)) reduce to 

r = ±E, (17) 



i = -(r 2 +a 2 + l 2 )^, 
aE 



(18) 
(19) 



The radial coordinate is uniformly described with re- 
spect to the affinc parameter while the equations gov- 
erning t and ip> are 



dt r 2 + a 2 + l 2 

~dr = ± A : 

dp ^a 

dr ~ A' 



(20) 
(21) 
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3 Energy extraction by Penrose process for 
Kerr-Taub-NUT spacetime 

Let us continue our assumption that the deflection in 8 
direction is reasonably small and orbits of the particles 
are in the quasi-equatorial plane 9 = it /2 + 86{t). Using 
the assumptions mentioned in previous section one can 
easily rewrite the expressions ([5]) and ([6]) in the approx- 
imation 0(89 2 (t)) as a quadratic equation in energy: 



aE 2 -2[3E + 7 + -{p r f + S(/) 2 + m 2 = 0, 

where we have used the following notations 
Mr + 1 2 ' 



4a- 



1 - 2 
Mr 



r 2 + 1 2 



r 2 + I 2 

2 

L, 
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a 2 + I 2 + 2a 



,Mr + l 2 
r 2 + l 2 



(22) 

(23) 
(24) 
(25) 



From the equation (f2"2"j) one can easily obtain the equa- 
tion of radial motion in the following form: 



r 2 = E 2 



(26) 



and the notation 



E 2 - 2 
E 2 (r 2 - 



(Mr + 1 2 ){L - aE) 2 

£(r 2 +/ 2 ) 
f a 2 + I 2 ) L 2 + A 



(27) 



denotes the effective potential of the radial motion of 
the test particle around rotating compact object with 
nonvanishing NUT parameter. 

In the Fig. [5] the radial dependence of the effec- 
tive potential of radial motion of the massive test par- 
ticle has been shown for the different values of the di- 
mensionless parameter Here for the energy and mo- 
menta of the particle the following values are taken: 
E/m = 0.9, L/mM = 4.3. The presence of the param- 
eter / slightly shifts the shape of the effective potential 
down. 

Now, as the particle falls through the horizon, the 
mass of the black hole will change by SM = E. There 
is no upper limit for change of mass of the central black 
hole. The infalling big number of particles with positive 
energy can essentially increase the mass of the black 
hole. But there is a lower limit on SM which could be 
added to the black hole corresponding to m = 0, p e = 
and p r = 0. Evaluating all of the required quantities at 
the horizon r — r + , we get the limit for the change in 
black hole mass as 



SM = 
where 
D = 4 



L z a(Mr+ + 1 2 ) 
r\ - 2Mr+ 



D. 



(28) 



1-1 



1 - 



a 2 (r\ +2Mr+ +31 2 ) 
(r\ + l 2 ) 2 



(r 2 + - 2Mr+ - I 2 ) (Mr+ + l 2x 



16a 2 



4 + 1 2 



(29) 



To be able to extract energy from the black hole 
(SM < 0), wc must therefore have 



L z a(Mr + + I 2 ) 
rl - 2Mr+ 



D > 0. 



(30) 



In the Fig. [3] the dependence of the extracted energy 
from the black hole on the small dimensionless param- 
eter I has been shown: with increasing the parameter 
I the relative extraction of the energy becomes more 
stronger. The graph shows that the extraction of the 
energy is directly proportional to the parameter I , that 
is with increasing Z, the extracted energy also increases. 



4 Conclusion 

We have studied the properties of the ergosphcre of 
the black hole in the Kerr-Taub-NUT spacetime. The 
dependence of the shape of the ergosphere from small 
dimensionless NUT parameter shows that the radius of 
the event horizon becomes larger. However the relative 
volume of the ergosphcre is decreased and the extracted 
energy from the black hole is raised up by the small 
dimensionless parameter I . 
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Fig. 2 The radial dependence of the effective potential of 
radial motion of the particle for the different values of the 
dimensionless NUT parameter I . 




Fig. 3 The dependence of the extracted energy from the 
black hole on the dimensionless NUT parameter I . 
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